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Abstract 

We investigate the stability of self-gravitating spherically symmetric anisotropic 
spheres under radial perturbations. We consider both the Newtonian and the 
'. full general-relativistic perturbation treatment. In the general- relativistic 

^ I case, we extend the variational formalism for spheres with isotropic pressure 

O i developed by Chandrasekhar. We find that, in general, when the tangential 

2 ' pressure is greater than the radial pressure, the stability of the anisotropic 

■ sphere is enhanced when compared to isotropic configurations. In particular, 

Q^. anisotropic spheres are found to be stable for smaller values of the adiabatic 

^ ', index 7. 
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H ■ 1 Introduction 



In a recent paper [1], we presented a series of new exact solutions of the 
Einstein field equations for self-gravitating, general-relativistic spheres with 
anisotropic pressure. We have found that the presence of pressure anisotropy, 
for a large variety of ansatze for its functional form, has important implica- 
tions for the physical properties of self-gravitating objects. Namely, both 
the object's critical mass and surface redshift are modified, and may violate 
well-known bounds for isotropic objects {2M/R < 8/9 and Zg < 2). We have 
shown that this is true not only for stars of constant energy density, but also 
for objects with p oc often used to model neutron star interiors. Given 
the fact that pressure isotropy is an assumption not required by imposing 
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spherical symmetry, it is clearly of great relevance to investigate if, indeed, 
these anisotropic configurations are stable against radial perturbation and, 
thus, better candidates to exist in Nature. 

The aim of the present paper is, then, to develop a formalism which 
can be used to test the stability of anisotropic spheres against small radial 
perturbations. Our formalism is a generalization of the variational principle 
used for investigating the stability properties of isotropic objects. We reduce 
the stability analysis to an eigenvalue problem, where the eigenvalues are the 
frequencies of the radial modes. 

The dynamical stability of isotropic spheres has been extensively studied 
by various authors [2, 3, 4, 5, 6]. A calculation based on the concept of 
extremal energy was presented by Fowler [3]. Cocke [4] performed a calcu- 
lation based on the method of extremal entropy. Gleiser, and Gleiser and 
Watkins applied Chandrasekhar's variational method to investigate the sta- 
bility of boson stars, self-gravitating spheres of complex scalar fields, which 
are naturally anisotropic [6]. 

Chandrasekhar considered the dynamical stability of isotropic spheres as 
an eigenvalue problem [2]. He used an analytical approach to compute the 
eigenfrequencies of radial oscillations for isotropic spherical stars. The study 
of the stability of a star thus becomes a Sturm-Liouville problem. 

The main result of these studies is that, for dynamical stability in general 
relativity, isotropic spheres must have an adiabatic index (or exponent) 



where k is a number of order unity, that depends on the structure of the star, 
and M and R are the star's mass and radius, respectively. For white dwarfs. 



The stability of anisotropic spheres in general relativity was studied nu- 
merically by Hillebrandt and Steinmetz [7]. An analytical approach in the 
spirit of Chandrasekhar's work for isotropic spheres, however, docs not seem 
to exist for anisotropic spheres. Our goal is to obtain this approach. 

This paper is organized as follows. In the next section, we obtain exact 
solutions for several examples of Newtonian anisotropic spheres, and study 
their stability properties. This will give us some insight into the effects of 
anisotropy on the stability of self-gravitating objects. We then proceed, in 
section 3, to derive the full general- relativistic perturbation formalism for 
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anisotropic spheres. In section 4, we apply the formahsm to anisotropic 
spheres of constant energy density. In section 5, we apply it to anisotropic 
spheres with p oc In both sections, we follow the exact solutions derived 
in [1]. In section 6 we summarize our results, and discuss possible avenues 
for future work. 



2 Newtonian Anisotropic Spheres 

2.1 Exact solutions for Newtonian Anisotropic Spheres 

We consider the dynamics of anisotropic spheres under the influence of New- 
tonian gravity. The equation of hydrostatic equilibrium with anisotropic 
pressure in Newtonian gravity is 

, m(r)p(r) 2 , , 

Pr--^^ + -iPt-Pr) (2) 

where is the radial pressure, pt is the tangential pressure, p is the energy 
density and 

m(r) = 47r / p{r')r'^dr', (3) 





is the mass contained in a sphere of radius r. This equation may be obtained 
as the Newtonian limit of the generalized Tolman-Oppenheimer-Volkov equa- 
tion for general relativistic hydro dynamical equilibrium, or it may be derived 
using the principles of Newtonian fluid mechanics. 
The pressure in isotropic spheres with constant density, Po, is given by 

Pr=^-^Po'{R'-r'). (4) 

We note that, in Newtonian gravity, the pressure at the center of a sphere 
with constant density and isotropic pressure can only become infinite if the 
radius of the sphere is infinite. 

We will now solve eq. 2 for various ansatze connecting pr and pt at con- 
stant density po- These ansatze are chosen so as to correspond to the choices 
we will make for the full general relativistic cases. 

Case I: pt — pr = C pi 

This ansatz assumes that the anisotropy term in eq. (2) is proportional to 
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the first term on the right hand side of the equation, i.e., the anisotropy is 
chosen to mimic the behavior of the purely gravitational term. This ansatz 
can be interpreted as the Newtonian limit of the ansatz that Bowers and 
Liang used to solve the full general relativistic TOV equation [12]. With this 
ansatz eq. (2) becomes 

p'r--Ypl^ + 2Cplr (5) 



and the solution is 



Since we are considering spheres with constant energy density po, from eq. 
(3), 
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m{r) = -npor^ . (7) 
Therefore, here we can also write 

Comparing this solution with the isotropic case {C = 0), we see that the 
radial pressure has the same spatial behavior in both cases, and they differ 
only by a multiplicative factor that depends on the amount of anisotropy C. 
If we define an effective density: 

then we can write 

Pr = ^p'(R'-r'). (10) 
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Thus, in this model the effect of the anisotropy can be considered as a scaling 
of the density of the sphere. We can take this scaling interpretation a step 
further by reintroducing into the expression for pr the gravitational constant 
G: 

Pr-plG[^-c){R'-r'). (11) 

We now define 

G^{i-'^)g (12) 
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Figure 1: Radial pressure Pr as a function of r for the ansatz Pt—Pr = Cp^r'^, 
parameterized by values of C. 

and find that 



Thus, we can also interpret anisotropy as a variation of the gravitational 
constant. This explains why positive values of C and hence smaller values of 

G lead to smaller values of the radial pressure. For smaller G the gravitational 
force between the particles in the sphere is decreased and this leads to a 
decrease in the radial pressure needed to stabilize the sphere. Negative values 
of C have the opposite effect, i.e., G is increased and correspondingly p^ is 
also increased. We plot the radial pressure Pr as a function of the radius r 
for several values of the anisotropy G in figure 1. 

We note that, for C = 27r/3 ((5 = 0), the radial pressure vanishes and 
becomes negative if C > 2n/3; in this case, no bound solutions are possible. 
It is interesting to note that the solution with G — 27r/3 has the following 
form 



Hence, for this particular solution, the sphere is held together by purely tan- 




(13) 




(14) 
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gential stresses. The particles that constitute the sphere are considered to 
be in random circular orbits [11]. 

Case II: pt — Pr — C po Pr 



The solution to eq. (2) with this ansatz is 



2tt 

Pr = ^PO 



For this solution, the pressure at the center is 



2tt 

Pc = ^PO 



1 - e 



and all values of C are allowed. For C small, 

2n 



and 



Pt = Ypl(R"-^")(^ + CPor'). 



Case III: Pt—Pr^Cp^ 



(15) 

(16) 

(17) 
(18) 



Here the solution has two distinct forms depending on whether C is 
greater than or less than zero. For C < 0, the solution is 



with 




Thus Pc becomes infinite \i E? = ^ (^) ' ^ ^^^sult that is not possible for 
Newtonian isotropic spheres with constant density. However, we note that 
the values of C for which pc becomes infinite in Newtonian gravity are quite 
large, on the order of 10'^ for 2M / R = 0.05. We will see that for the general 
relativistic case when 2M/R ~ 1, the values of C for which p^ becomes 
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Figure 2: Core pressure Pc ^ a function of C for the ansatz Pt—Pr = Cp^'^r'^, 
parameterized by values of 2M/ R. 

infinite are of order 1. We plot the core pressure Pc as a function of C for 

various values of 2M/R in figure 2. 
When C > the solution is 



27r\V2 



exp 




) Po [r' 




- 1 


exp 




1 Po 




+ 1 



(21) 



For this solution pc is always positive and finite. 

2.2 Stability of Newtonian Anisotropic Spheres 

We now proceed to investigate the effects of small perturbations on the solu- 
tions obtained above. This study is important since it allows us to compute 
the frequencies and normal modes of oscillations, enabling us to establish 
the dynamical stability of our solutions. We will follow closely the formal- 
ism outlined in Shapiro and Teukolsky [8] for isotropic spheres, modifying it 
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where necessary to the anisotropic case. 

It is useful to distinguish between Eulerian and Lagrangian perturbations 
of the fluid variables. If Q{x, t) is any fluid variable, the Eulerian change 
relative to the unperturbed value Qo{x, t) is defined as 

dQ{x,t) =Q{x,t) -Qo{x,t). (22) 

The Lagrangian change is defined as 

AQ{x, t) = Q[x + i{x, t),t]- Qo{x, t) (23) 

— * 

where ^(5;, t) is an infinitesimal displacement of the fluid element. 

The Eulerian approach considers changes in the fluid variables at a par- 
ticular point in space, whereas the Lagrangian approach considers changes 
in a particular fluid element. The relationship between the two is 

A = (5 + f • V. (24) 

The following equations govern the dynamics of the unperturbed system: 

1. The continuity equation that connects the density p and velocity v, 

% + {pv)' = ; (25) 

2. The momentum equation 

dv 1 , ^,2, , 

-f,=—Pr-^ +-iPt-Pr) , 26 

at p r 

where 

d d d 

3. Poisson's equation, the equation that determines the gravitational po- 
tential 



-(r2$')' = 47rp . (28) 
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We have adopted spherical symmetry since we are considering radial pertur- 
bations. 

A Lagrangian perturbation of the momentum equation gives, 



We note the following: 

• Prom the continuity equation it follows that 

• A perturbation of Poisson's equation gives 



(29) 



(30) 



(5$)' = -47rpe (31) 
• The adiabatic exponent 7 is defined using the following expression, 



p 



(32) 



• We also find it convenient to introduce the following symbol, 

Xl = Pt-pr. (33) 

Using eqs. (27)- (30) we can evaluate each term of the perturbed momentum 
equation: 



. dv d^^ 
A — = — 
dt df^ 



(34) 



/ 1 A / 2 ^ , 1 
-—Pr + -^Pr = —ZPr + ' 

p pr p 



(35) 



2 d$ 2e , 2 
r dr pr pr^ 



(36) 
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and 

A fin) = -4(Cr^)'n - + Am. (37) 

Combining all terms in eqs. (34)- (37), we find that radial perturbations are 
governed by the following equation: 



P 



4 6 2 2 

-CK-— ^n--(e)'n--cAn 

pr pr"^ pr pr 



(38) 



We now assume that all variables have a time dependence of the form 

e*"^*. Substituting this form of the time dependence in the above equation, 
we arrive at an eigenvalue equation for radial oscillations of a Newtonian 
spherical star; 



wr-,{r'0'\ - -U + + -(O'n + -CAn + pu'^ = o . (39) 

The boundary conditions for this equation are 

^ = at r = 0, (40) 

Ap^ = at r = R. (41) 

Equation (39) subject to the boundary conditions (40) and (41) is a Sturm- 
Louville eigenvalue problem for a;^. The general theory of these equations 
gives the following results [10]: 

1. The eigenvalues are real and form an infinite discrete sequence. 



2. The $,n are orthogonal with a weight function pr^: 

inimpr'^dr = 0, m^n. 

3. The form a complete basis for any function satisfying the boundary 
conditions 40 and 41. 
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An important consequence of these results is that, if the fundamental 
mode of the star is stable {uq > 0), then all radial modes are stable. Con- 
versely, if the star is radially unstable, the fastest growing instability will be 
via the fundamental mode {u!q more negative than all other cu^). 

Equation (39) can be solved for cu'^. Multiplying by ^r^ and integrating 
from to i? we find 

/J'{7P.^W + +4rey-6fn-2ree'n-2reAn}dr ^ ^ 
^ = 5 — . (42) 

We will now compute the frequency of oscillation for some anisotropic 
spheres. We will only consider models with 7 = const. First, we note that 
for stars with isotropic pressure and constant density, under a self-similar 
deformation ^ = const x r, 

^2 = 47rpo(7 - I) ■ (43) 

This result says that isotropic stars with 7 = 4/3 are marginally stable. If 
7 is less than 4/3 then dynamical instability will occur, while if 7 is greater 
then 4/3 the star is stable relative to the deformation = r. 

Computing uj from eq. (42) for the anisotropic model Case I, with the 
deformation ^ = r, we find 



LO^ = 6po 



(44) 



Thus, not surprisingly for this model, anisotropy results in a scaling of 
the frequency of oscillation. We had already seen earher that the effect of the 
anisotropy was equivalent to a scaling of the density, and since the frequency 
of oscillation is proportional to the density, this result is expected. Thus, for 
this case, positive anisotropy may slow down the growth of instabilities, but 
will not reverse their trend (recall that we must have C < 27r/3.) 

We next consider the model Case II. For this model. 



4 16^ 

7-5 + 5jCp„fl^ 



(45) 



The fundamental frequency occurs for a;^ = 0, and this corresponds here to 
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Figure 3: Adiabatic index 7c, for which uP' = 0, as a function of C for the 
ansatz Pt — Pr — CpoPrV^ (Case II) parameterized by values of p^R. 

7c = ^ - f^CpoR' . (46) 

We observe that, since poR^ is always a positive quantity, depending on 
the sign of C the value of 7 for which u;'^ — can be less than or greater 
than 4/3, indicating that positive (negative) values of anisotropy can increase 
(decrease) the stability of the star. 

We plot 7c as a function of C for u;^ = in figure 3 for the ansatz 

Pt-Pr = CpPrV'^. 

This concludes our stability study of Newtonian anisotropic spheres. We 
have seen that the presence of anisotropic pressure in a self-gravitating sys- 
tem can have dramatic effects on the dynamics and stability of the system. In 
particular, there are some novel features that are present only if the pressure 
is anisotropic, e.g., infinite core pressure, zero radial pressure and stable ob- 
jects with 7 < 4/3. We will now proceed to study the perturbation problem 
for relativistic anisotropic compact spheres. 
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3 Stability of General Relativistic Anisotropic 
Spheres: General Formalism 

In this section we will study perturbations of exact solutions of the gen- 
eral relativistic field equations for anisotropic spheres. In particular, we will 
be concerned with perturbations that preserve spherical symmetry. Under 
these perturbations, radial motions will ensue. We will develop an analytical 
approach generalizing work by Chandrasekhar for isotropic spheres [2]. 

In this section po, z/q, Xo,Pto Pro values of the dynamical variables 
that satisfy the equations for static equilibrium. The perturbed variables will 
be written as p, i^, A,pt and Pr, respectively. 

3.1 The Perturbed Energy-Momentum Tensor 

The energy-momentum tensor for a spherically symmetric spacetime is [9] 

={P + Pr)uX - 9> - 1^.1" {.Pt - Pr) - Kk^'ipt - Pr) . (47) 

Here, 

'"' = ^- («' 

and 

1,^5% r^5l k,^5% k^^St (49) 
Since we are considering only radial motions, we will take 

= e"^, Ut = e^, (50) 

and 

li = 2 ^ ve ° 2 , (51) 

with 

-I- 

It should be quite clear when a subscript refers to a time-like coordinate, 
such as li*, or tangential pressure, such as pt- We reserve the index o to 
unperturbed metric and physical quantities. 

Writing p = po + Sp, Pr = Pro + (^Pr, Pt = Pto + Spt, A = Ao + SX, 

and u = Vo + we find that, to first order in v 
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T} = p, T; = -pr, Tl = Ti = -pt (53) 

and 

T[ = iPo + Pro)v, Tl = {p, + pJve^>^''-^^\ (54) 

3.2 Perturbations of the Dynamical Variables 

The set of equations governing radial motions can be written as [2]: 

(re"^)' = 1 -87rpr% (55) 
I/' = r(e^ - 1) + Sirprre^ , (56) 

^e-^ = SttT; , (57) 
fl + (X + u)^ -p/ - 1 (p + p,) i.' + (58) 

with 

U={pt-pr). (59) 
The zeroth order (or static equihbrium) equations are: 

(re-^°y = 1 - Snpor'^ , (60) 

= r(e^° - 1) + Sirprore^" , (61) 
and ^ 

Pr'o^-{Po+Pro)^ + l^o. (62) 

We also have the identity 

{K + <) = 87r{pro + Po). (63) 



Wc now linearize eqs. (55)- (58), taking into consideration eqs. (60) - 
(62). Since we consider all perturbations to be of order v, we find that to 
first order in v, eqs. (55)- (58) imply 
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(5u' - vj5\) = — (5A + 87r% , 

(5A^ = -Stt (po + Pro) 



(64) 
(65) 
(66) 



(Pro + Po)v+ {Spr)' + T; iPvo + Po) i^^)' " o (<^P + ^Pr) K + = ^■ 

(67) 

(68) 



We now introduce a "Lagrangian displacement" ^ defined by 



dt' 



Integrating eq. (66), we find that 



SX- 



-^-87r{po + Pro)C- (69) 
Taking into consideration eq. (63), the above equation becomes 

5X = -e (A; + u',) . (70) 
We can now combine eqs. (64) and (70) to get 

Sp--^lr'(Po+Pro)^]' ■ (71) 

Substituting for the expression for p^o from eq. (62) into the above equation, 
we find that 



Sp = - iPo+Pro) ;^e-/^ (r^e-^/^e) - f . 

We now consider eq. (65). Using eqs. (62) and (63) we find 

1 



(72) 



{Pro + Po) (H' 



SPr - {Pro + po) {K + - ) ^ 



{\o + iyo) . (73) 



We note that eqs. (70), (72) and (73) allows us to express 6X, dp, and 
Su in terms of 6pr, v, and the unperturbed variables. We need to impose an 
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extra condition on the system in order to obtain an expression for 6pr. The 
condition we shall impose is the conservation of baryon number. Further, we 
note that 511 can always be expressed in terms of the unperturbed variables 
once Spr is given in terms of these variables. 

Chandrasekhar [2] derived an expression for 6pr from the law of conser- 
vation of baryon number in general relativity. Since we arc not making any 
new assumptions with respect to Chandrasekhar (except, of course, we are 
considering anisotropic pressure) , we will only outline the basic steps of his 
derivation here. 

The law of conservation of baryon number density in general relativity 
can be written as 

(nM");a = , (74) 

where n is the number density for baryons and is the four-velocity of the 
fluid. Taking 

n = no{r) + 5n{r, t) , (75) 

and recalling that, to first order in u"" is given by eqs. (50) and (51), eq. 
(74) becomes 

e-''°'^5n + 4 (nor^ve"^"/^) + -noe""''^5\ + -e-^'-'^v (A^ + !/„)' = 0. (76) 
Since v — ^, eq. (76) integrates to give 

5n + ^ {uor'^ve-'"''^)' + [5X + ^ (Ao + Vo)'] = 0. (77) 

The last term on the left-hand side of eq. (77) vanishes on account of eq. 
(70), and we obtain 

5n = (norV-"°/2j . (78) 

The first law of thermodynamics in general relativity is obtained by combin- 
ing 

u.T'^'',^ = 0, (79) 

with the law of conservation of baryon number given by eq. (74). Using the 
expressions for T^" and from eqs. (47), (50) and (51), we find that, 

Prd(-\+d{^\ + -{pt-pr)-dr^Q. (80) 
\nj \nj r n 
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Thus, in general, the equation of state is given by 

n = n{p,pr,'n.) . 



(81) 



However, since we are considering systems where the tangential pressure is 
given in terms of the radial pressure and the density (recall that in generating 
exact solutions in the previous chapter we assumed various ansatze for pt — 
Pr), we will take 

n = n{p,pr) . 

For this n, we have 



bn = ^— dp + ^— dPr + -^dr. 
op opr or 

Substituting for 5n from eq. (78) and 5p from eq. (72), we find that 



(82) 
(83) 



,dno 



(84) 



ar op r op opr or 



Dividing through out by (drio/dpr) gives, to first order, 

- i-o5 (rV"/^0' = -(P.. + Po)^(rV"/^04 (85) 



dr dno'"" J.2 

OPr 



dpr 



1 dup „^ dno 

Opr Opr Opr 



Solving for 5pr we find 



5pr = -Pr'i 



1 



Pro dp^ 

We can rewrite this as 



dn, 



- -Q^iPo+Pro) 



Spr = -PrJ - Wro^ {r'e^"^'^)' + 

with 7 being the adiabatic exponent defined as 



2^^ dpro 

r" ^ ^ r dpo ' 



r dpo 
(86) 

(87) 



7 = 



1 



Pr{dn/dpr 



n - {p + Pr) 



dn 

dpr 



(88) 
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3.3 The Pulsation Equation 



We now assume that all perturbations have a time dependence of the form 
^iu)t_ Further, considering 5A, 5z/, 5p, 5pr and 511 to now represent the am- 
plitude of the various perturbations with the same time dependence we find, 
from eq. (67), that 



(89) 



Substituting the expressions for the various amplitudes in the above equa- 
tion we find 

{Po + Pro) ie^"""" = -iProi)' -\iPo+ Pro) + ^) C (A^ + l^'o) (90) 



_g-(Ao+2i/o)/2 



g(A„+2.„)/2^^^^^^^2g-.„/2^y 



+ [l^K + K]PrJi 



,(Ao+2i/o)/2 2^j-^^'9pr 

r "dp 



r 



Prom eq. (62) it follows that 



,2, 



= -(p; + (Po + (^n,)' , 



and 



-2rp/„ + 4n„ 



Also, we have 



(91) 
(92) 

(93) 



Using eqs. (91), (92) and (93) in (90) we arrive at the pulsation equation 
i.e., the equation that governs radial oscillations: 
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(94) 



87re^''(no+p,J(p„+pje- 



{po + Pt 



{ProTi + 



APo^Pt^ r^{po+Pro) 



-(A<,+2i/o)/2 



,(A„+2.„)/22^n, f ^ + 1' 



The boundary conditions imposed on this equation are 

^ = at r = and Spr — a,t r — R. 



(95) 



The pulsation eq. (94), together with the boundary conditions eq. (95), 
reduce to an eigenvalue problem for the frequency cu and amplitude ^. This 
is equivalent to the Sturm- Liouville problem we encountered while studying 
Newtonian gravity in chapter 2. Multiplying eq. (94) by r'^^e^^'^'^y^ and 
integrating over the entire range of r we find that 



Jo Jo 



(96) 



Jo 

rR 



(A+3i/)/2 



R 

- I e 



^2f2 

(A+i/)/2 ' S ^„/^2 



iPrYdr 



Stt / e^^^+''y^(U+pr)(p + PryC^dr + 4 f e^^+'' 
Jo Jo 

~p 



(P + Pr) 



)/2. 



(p + Pr) 



Jo ip + Pr) Jo Jo 



Jo 



'u/2 



g(A+2.)/2 2^n f^ + l 

r \dp ^ 



r^^dr , 



where we have dropped the subscripts as no longer necessary. The orthogo- 
nality condition is now 



R 



,(3A-i/)/2 



(p + Pr)rW^^r = (i^i) , 



(97) 



where and are the proper solutions belonging to different eigenvalues 
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4 Stability of Anisotropic Spheres: p = const. 

Chandrasekhar studied the dynamical stability of isotropic spheres with con- 
stant density p and constant adiabatic index 7 using the above formalism. 
The full solution of Einstein's field equations for constant-density isotropic 
spheres is well known [9]. Thus, writing 



= l-rf 



and yi = 1 - ^ = 1 - 



with 



87rp' 



the complete interior static isotropic solution for p = const is 

y-yi 



p 



3z/i - 



\ and e" = ]-{3yi - yf 
y 4 



(98) 
(99) 

(100) 



Here, we will apply the formalism just developed to the Bowers-Liang solution 
for anisotropic spheres [12] , 



Pr 



y'^ - yl"" 
2>y? - y'^ 



1 and e^^-i?,yl^-y'Qf''^ (101) 



and 



with 



^y'^yi 



2Q 



y^ {Zyl^^ - y^Qf' 



1 3C _ 1 k 



(102) 



(103) 



In his stability analysis for isotropic stars, Chandrasekhkar used the following 
trial function 

^ = rje''/^ = lrj{3y^-y) , (104) 
and found that, to first order in 2M/R, the frequency of oscillation is given 

by 

^' = :r^[(37-4)--^-(^)(547-53)]. 



2o;2' 



14' R 



(105) 
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The first term reproduces the results for Newtonian isotropic spheres [ cf. eq. 
(43)] and the second term represents a correction due to general relativity. 

We now turn our attention to the anisotropic case. A trial function that 
generahzes (104) to include the effects of anisotropy is 



rje 



Qv 



(106) 



However, we found that for the anisotropic case the corresponding inte- 
grals in the expression for uo"^ cannot be computed analytically if this substi- 
tution is made. In order to compute the integrals analytically, we used trial 
functions of the following form 



6 = V'^e'i^ , 



(107) 



and 



(108) 

For the trial function ^i, we found, after integrating all terms in eq. (96), 
[for small 2M/R] 
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Thus, for this model, stable oscillations will occur if 
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For the trial function ^2 integrating eq. (96) gives. 



49 
32^ 



245 



128^2 



48 931 586 
.^~49~ ^490^ "245! 



2M 
(109) 



(110) 



2M 



and the condition for stable oscillations becomes 
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An examination of the two expressions for 7 above shows that for positive 
anisotropy {k > 0) 7 is smaller than the corresponding isotropic value, 
implying that positive k leads to more stable configurations, while negative 
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values of k will have a destabilizing effect. Further, we note that since for 
small values of k the various new analytical solutions we found previously for 
constant-density anisotropic spheres [1] have a similar form to the Bowers- 
Liang solution, we expect the relationship found here between the sign of k 
and the stability of the sphere to hold also for those new solutions. 

It is useful to compare the results Chandrasekhar found using the trial 
function ^ [eq. (104)] (denoted by jch) , with the results obtained with our 
trial functions and ^2 (denoted by 71 and 72) in the isotropic limit {k = 0). 
For stable oscillations we must have 

^ |(HM).,333 + 0.4623(H^), (113) 
3«'2^''= 1.280 + oW' 
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3 + 






32 


71 
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25 ' 






64 


72 


> 


49 " 



75' i? ' ■ ' i? ' ' 

20, 2M, ,2M, 
_(_) = 1.306 + 0.4081(-). 

It is known that Chandrasekhar's trial function becomes exact in the limit 
{2M/R) — > 0. Comparing the results for the three trial functions above, we 
see that in the exact limit our results differ from the exact value by ^ 5%. 
This leads us to believe that our trial functions generate results that are 
qualitatively correct. 



5 Stability of Anisotropic Spheres: p oc 

In [1] we found several exact solutions for anisotropic stellar configurations 
with the following expression for the energy density 

P=8^(^ + »). (114) 

where both a and b are constants. The choice of the values for a and b is 
dictated by the physical configuration under consideration. For example, 
a — 3/7 corresponds to the Misner-Zapolsky solution for ultra high-density 

neutron star cores [13]. 

If we model the pressure anisotropy as 

^*"^^^^(;^ + ^) ' ^^^^^ 
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then an exact solution of the field equations with b = d = is 

e-^^l-a = 7o (116) 

with 

.-^ii^. (118) 

(1 -a)2 

and the constants and A_ fixed by boundary conditions. For the case 
under consideration here {b — d — 0), the boundary conditions are 

^-m ^ ^.(R) ^ p e^(^)^|« = - . (119) 

dr R 

Applying the boundary conditions we find 

= | + and A.^A4q^-q). (120) 

The radial pressure for q real, after substituting the expressions for A+ and 
A_, is 



r2 [ (3/0^ - 1 + 2qll)B?<i + (1 - 3/0^ + 2qll)r'^^ 

(121) 

For this case we found that using the following trial function 

e = r2(p + p,)e'^ (122) 

all the integrals were, after some tedious work, exactly integrable. In table 

I, we present results for the frequencies of radial oscillations cj^ as a function 
of the anisotropy parameter, c, for given values of the density parameter. 
We also give, in table II, the values of 7c above which stable oscillations are 
possible. Here we see that the effect of a positive anisotropy is to reduce the 
value of 7, thus giving rise to a more stable configuration when compared with 
the corresponding isotropic model. In particular, for the Misner-Zapolsky 
solution (a = 3/7), we find that a small positive pressure anisotropy in the 
equation of state improves the neutron star's core stability. 
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a = 2/9 


u'^E? = 0.95(7 -1-79) + (101.1 -52.6 7 )c 


a = 2/7 


cj2/?2 = 2.3(7 - 1-83) + (122.3 - 59.3 7)0 


a = 3/7 


uj'^R^= 0.57(7 - 1-93) + (15.2 - 5.l7)c 


a = 3.4/7 


uj'^K^ = 0.4(7 - 2.6) +(8.9 - 2.3 7)0 


a = 3.49/7 


uj^R'^ = 0.36(7 - 2.76) + (8.0 -1.97 7)0 



Table 1: o;^ vs. c for given values of a. 



a = 2/9 


Cmax = 0.0016 


7c = 1.79 -6.87 c 


a = 2/7 


Cmax = 0.0028 


7c = 1.83 - 13.39 c 


a = 3/7 


Cmax = 0.083 


7c = 1.93 - 5.55 c 


a = 3.4/7 


C-max 0.11 


7c = 2.6 - 2.84 c 


a = 3.47/7 


C-max 0.12 


7c = 2.75 - 7.29 c 



Table 2: 7c c for given values of a. 



6 Conclusion 

We have studied the stability of anisotropic compact spheres against radial 
perturbations in the framework of Newtonian gravity and general relativity. 
In both cases we have seen that the presence of anisotropic pressure can have 
significant effects. 

We have found that there arc Newtonian anisotropic spheres with con- 
stant energy density whose core pressure can go to infinity, without requiring 
the radius of the sphere to be infinite. A result of this nature is not possible 
for Newtonian isotropic spheres with constant energy density. Furthermore, 
a stability analysis of some of these models shows that there can exist stable 
anisotropic spheres with an adiabatic exponent 7 < 4/3. In the correspond- 
ing isotropic case, instability immediately sets in if 7 < 4/3. 

We have extended the formalism developed by Chandrasekhar to study 
the stability of general relativistic isotropic spheres against radial perturba- 
tions to anisotropic spheres. In particular, we have applied this formalism 
to study anisotropic spheres with constant energy density and with energy 
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densities with an profile, used to model ultra-dense neutron star inte- 
riors, for example. We have found that in both cases there can exist stable 
relativistic anisotropic spheres with values of the adiabatic exponent that 
would necessarily imply instability in isotropic spheres. In particular, this 
is true whenever the tangential pressure is larger than the radial pressure 
for all models we investigated. These results may explain the higher stabil- 
ity of certain neutron stars with anisotropic deviations near their core, and 
other gravitationally-bound compact objects such as boson stars, which are 
naturally anisotropic. Work along these lines is currently in progress. 
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